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Abstract 

We have studied the energy and persistent current spectra for a relativistic fermion on a 
one dimensional ring in different mass regime. We have thoroughly analysed the nonlinear 
nature of the persistent current and its dependence on particle mass and ring radius. For 
a particular radius of the ring we find the existence of a critical mass at which the current 
becomes independent of flux, the value being same as the current produced by a massless 
particle. This value becomes the lower cutoff of the current as we further decrease the mass. 
The total current spectrum on the other hand shows periodic plateaus and becomes discrete 
as we go below the critical mass. Finally we study the susceptibility of the system and 
showed in details how it changes its signature while going from a higher mass to lower mass 
limit. 

1 Introduction 

The existence of a persistent current in a normal metallic ring threaded by a magnetic flux 
has stimulated a large number of studies in last few decades pQ-[9]. The current is periodic 
in flux like the current in a superconducting ring |10|-[13j. Due to the recent advances in 
fabrication techniques, it has become possible to observe this current experimentally |14j-|17j 
in mesoscopic systems. Such current is also found in carbon nanotube rings [18j-|20j. In 
these carbon based structures electrons behaves as massless Dirac particles and they require a 
relativistic description. Although there are a good number of fine literatures on the energy levels 
and persistent current of relativistic fermion regarding graphene |21j-[26j. topological insulator 
|27| and for Dirac electron [2H] , there is still lack of a complete description of persistent current 
for relativistic fermions and its variation with particle mass. Besides, it is well known that 
superconducting as well as mesoscopic rings are perfect diamagnets, but carbon nanotube tori 
[29j are found to posses a positive susceptibility. There is no vivid discussion on how the 
transition occurs. 

In this paper we have studied the nature of energy and persistent current spectra for a 
relativistic particle on a one dimensional ring. Our scheme is valid for all the mass regime and 
can explain the massless case as well. For heavy fermions our finding matches well with the 
nonrelativistic result [3j. Interesting features are revealed as we move towards the relativistic 
regime, ie. for low mass limit. The flux dependence of the current becomes more and more 
nonlinear as we decrease the mass. For a particular ring radius we find the existence of a critical 
mass at which the current becomes constant. The value of this current is same as the current 
produced by a massless particle and it appears to be the lower cutoff for current as we further 
decrease the mass. This nature is prominent for massless particle as well. We further study 
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the nature of the total current and find that it shows periodic plateaus. It becomes a multiple 
step function at critical mass and becomes discrete below that mass. The same periodicity 
is also observed in the susceptibility. We find that its amplitude change from negative to a 
positive value as we lower the mass which is further clarified by studying the behaviour of the 
susceptibility with respect to mass at a particular flux value. 

The organisation of the paper is as follows. In section 2 we give a brief description about the 
Dirac equation and its solution on a circle. Then we solve it in presence of a Aharonov-Bohm 
flux and obtain the energy spectrum. We derive the resulting persistent current in section 3 
and show its behaviour for different mass range specially for a massless fermion. In section 4 we 
discuss the behaviour of total current and in section 5, the nature of susceptibility is discussed. 
Section 6 contains the discussion and the proposed experiment to verify the result. 



2 Relativist ic particle on a one dimensional ring 

We start with a very brief description of Dirac equation on a circle and its solution. Then we 
will derive the solution in presence of an Aharonov-Bohm flux and obtain the energy spectrum. 
Dirac equation in two dimensions is given by |30] 

(c<7 • p + a^mc?)^ = Eip (1) 

where i(j t = (Ci? C2) is a two component spinor. a = a\x + a%y where a\, 02 and 03 are Pauli 
matrices. Due to the symmetry of the problem, it would be easier to work with polar coordinate 
[ip) in which a ■ p becomes 



a ■ p 




f>—i<p 9 




(2) 



where R is the radius of the ring. 

Now let us consider a perpendicular magnetic field piercing the ring (fig[[$). A suitable 
choice for vector potential (under symmetric gauge) is 



A 



2^R 



(3) 



where $ is the total flux passing through the ring. 




Figure 1: A one dimensional ring of radius R threaded by a magnetic flux <£. The field is finite 
only within the shaded region. 

Corresponding Dirac equation can easily be obtained by the replacement 

— — — 

d<p^ dtp +t <s> [ ) 

where $0 = — • The energy of the system is thus given by [28] 



^ = ^v+(^) 8 a+±)(i+±-i) (5) 
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Figure 2: Energy as a function of flux. Here we have used h = c = e = m e = 1 (e and m e 
being electronic charge and mass) in which R = 10 and m = 1. The energy minima are at 
$/$o = I + 1/2. This shift is a signature of relativistic correction. 



where I = 0, ±1, ±2, • • • . 

One can see (/zgG]) that the energy varies parabolically with respect to and has a peri- 
odicity of ^o- This is similar to the case of a nonrelativistic particle [3J. The only difference 
is the fact that the minima of the energy is at half integral multiple of <J>o- This shift by $o/2 
appears due to relativistic correction [28J . 



3 Persistent current for a Dirac particle 



We will now discuss about the nature of energy spectrum and resulting persistent current for 
different mass regime. 

Persistent current is the perpetual current flowing in a normal metallic ring in presence of 
a Aharonov-Bohm flux. It is given by the flux derivative of the energy [31 [13] 



I, 



(6) 



For a Dirac particle on a one dimensional ring whose energy is given by ([5]), the persistent 
current is 
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(7) 
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Figure 3: Variation of current with flux for R = 10 and m = 1. Here we have used the same 
unit system as figM 



One can regain the nonrelativistic result for m — > oo. At this limit the current is given by 

1 

2 



I, 



l \m— too 



eh f 
" 2-KmR 2 1 + W 



(8) 
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The expression looks similar to the nonrelativistic current [3]. The only difference is that 
like the energy spectrum, it also undergoes a shift by <3?o/2 with respect to the flux. However 
the behaviour of current at relativistic regime (m — > 0) is rather complicated. To have an idea 
about the flux dependence of the current we will study the susceptibility 0/(k, $) defined as 

«) = 7^ = ^~ ; — 77^2 ( 9 ) 
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where 



2-KmR 2 



h - 2 



mcR / 

Now we will discuss the behaviour of fiz(3>, k) for different range of $ and k. 



(10) 
(11) 



a. k — > : r2z(<3?, «) — > —Iq/^q which is negative, constant and independent of I. Hence the 
current will vary linearly with flux for all I values. This is actually the nonrelativistic 
(m — > oo) limit. 

b. k = 4 : S7;(^>, k) = 0, i.e. the current will not depend on the flux. We call the corresponding 
mass to be the critical mass for radius R, given by 

"* " in (12) 

The constant current is given by 

I sat = ±^ = ±— (13) 
k 2-kR v ' 

This is actually current generated by a massless particle with charge e. 

c. k > 4 : This is the low mass limit. At $ — > oo limit the current becomes constant, given 
by (USD- For k > 4, k) blows up at 



$ 



div 



which are actually the roots of the denominator of ([9]). The divergence in ri;(<J>, k) suggests 
a discontinuity in current. For k — > oo (m = 0, i.e. massless limit) the divergence occurs 
at 7?- = — I and £■ = — I + 1. 

The nature of energy and current in these three k region is shown in fig® Fig [2] and [3] show 
the energy and current spectra for k = 0.01 (m = 1) which is an example for n — > regime. 

One can also find the same result by keeping m fixed and varying R as a parameter. In that 
case there would be a critical radius given by 

R?ri = ^ (15) 

at which the current becomes constant. This is actually half of the Compton wavelength of 
the particle. However by tuning the radius it is not possible study the phenomena for massless 
particles. 
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Figure 4: Nature of the energy and current for different values of k. Here we use (H = c = e = 
m e = 1) R = 10 and vary the mass. We have taken m = 0.1 (k = 1), m = 0.05 (k = 4) and 
m = 0.01 (/c = 100). m = 0.05 is the critical mass (|12p at which energy becomes linear in flux 
and current attains a constant value I sa t (|13|) . For n > 4, J sa ^ is the lower cutoff of the current. 
This creates a gap of 2I sat is current spectrum for k > 4. 



3.1 Persistent current for a massless Dirac particle 

As a special case, we will discuss here the nature of persistent current caused by a massless 
particle. The energy of such particle can be obtained by putting m = in ([5]), and is given by 



Corresponding persistent current by definition (JS]) is 



ec 



1 ^ $0 2 



2vri? 



The variation of E? and I. with flux is shown in fig® 
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Figure 5: Nature of energy and current for massless particles. The lower cutoff value for current 
is given by ec/27rR. Here we used R=10 with the same unit as figWl- 
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From the denominator of (|17p one can readily see that for $/$o = ~~ I an d ^/^o = —I + 1 
the current diverges. This is something that does not happen for a nonrelativistic particle. At 
asymptotic region (^/^o ~~ * °°) the current becomes constant (fig\5§ and is given by (|13|) which 
is the current produced by a massless particle with charge e circling in a ring of radius R. 

4 Total current for odd and even number of rings 

Experimentally what one can observe is the total current for a large number of rings |14j-|17|. 
From (|17p one can understand that an exact analytic expression for total current is not possible 
for our cas^]. We numerically evaluated the total current for 1000 and 1001 isolated rings, each 
with radius R = 10 and for different m values. The result is shown in fig® 
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Figure 6: Total current for 1000 (Blue) and 1001 (Red dashed) isolated rings each with R=10 
(H= c = e = m e = 1). Mass and corresponding k value is written on the top of each plot. For 
large mass (a) total current depends linearly on flux. As we decrease the mass (6) the relation 
becomes nonlinear. For critical mass (|12fl . total current is a step function of flux (c), the step 
heights being multiple of I sa t (TTH|) . For lower mass (d), total current tends to diverge at values 
given by fjl4[> . the mean current being multiple of I sa t- On further decrease of mass (e) the 
spectrum becomes narrower and the divergences come closer to the integral flux values. As we 
go to the zero mass limit (/), total current appears as sharp pulses at integral values of flux. 

From /igj6]we see that for large mass, total current varies linearly with flux {figWp)- This 
is consistent with the nonrelativistic result. As we decrease the mass, relativistic effects starts 
dominating and the current starts showing nonlinear behaviour {fig®)). This is manifested 
by the periodic appearance of flat regions in figM). At m = rn^ ri (fT2|) (k = 4) total current 
becomes a multiple step function of flux (/i^Et), the step heights being integral multiple of 
Isat (fl3j) and each step is 2I sat high. As we further decrease the mass (k > 4) the total current 
starts showing divergences (fig^d). The point of divergence is given by (|14p . The average of 
the plateau between two divergences is at some multiple of I sa t an d located at integral values 
of flux. The width of the plateau decreases with the decrease of mass (/ifl'JSfe). As we go to the 
massless limit (fig^f), the total current appears as sharp pulses at integral values of flux. 

For an approximate result at large mass limit see |28| . 
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5 Susceptibility of the collection of rings 



The total susceptibility is evaluated by taking the flux derivative of the total current. It is well 
known that superconducting rings are perfect diamagnets and so are the mesoscopic rings. Let 
see the nature of the susceptibility for the relativistic fermions. 



(a) m = . 100, x = 1 . 
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Figure 7: Nature of susceptibility for a ring with radius R=10 at different mass which is written 
on the top of every plot. We have used the same unit as fig® For large mass susceptibility is 
negative (a) and almost constant. It shows periodic oscillation with a negative amplitude (b, c) 
for m > m cr i and becomes zero at critical mass ((f). For lower mass it again shows oscillation 
(e, /) but with positive amplitude. The blue and red dashed lines correspond even and odd 
number of electrons which overlaps on each other. 

The change of the magnetic susceptibility with mass is quite interesting for relativistic 
fermions (shown in figH}- For large mass the susceptibility is negative and almost constant 
{figWp)- The ring behaves as a diamagnet like a superconducting ring. As we decrease mass 
it starts showing periodic oscillation with negative amplitude (figHIb, c)) with period $o- The 
maxima are located at integral flux values. The susceptibility becomes zero at critical mass 
figWtf with a sharp divergence at half integral values of flux which are due to change of step 
height in flux-current relationship curve figWp- As we further decrease the mass, the periodic 
oscillation appears again, but this time with positive amplitude figWi^,f)- The susceptibility 
has a minimum value at integral flux and encounters a positive divergence at half integral 
values of flux. The positive magnetic susceptibility is also reported for carbon nanotube torus 
|29j where the electrons behave as massless Dirac fermions. To understand the the dependence 
of susceptibility on particle mass it is better to focus at the vicinity of an integral value of flux. 
Here we choose $ = and study the variation of susceptibility with respect to mass. 

The nature of change of susceptibility with respect to mass is clear from fig® For large mass 
we see that the susceptibility is negative and varies little with mass. This is the nonrelativistic 
limit where the rings behave as diamagnets. The susceptibility increases with the decrease 
of mass and becomes zero at critical mass. On further decrease of mass the susceptibility 
becomes positive and keeps increasing as we proceed to the zero mass limit. The collection 
of single particle fermionic rings thus makes a smooth diamagnetic to paramagnetic transition 
with decrease of mass. 
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Figure 8: Variation of susceptibility with respect to mass for a ring with R=10 at $ = 0. We 
have used the same unit system as figM 

6 Discussion 

In this paper we analyse the nature of energy and persistent current spectra of a relativistic 
fermion on a one dimensional ring. For a heavy fermion the result is similar to that of a 
nonrelativistic particle. The spectra start changing dramatically as we go to the low mass 
regime. For a given radius we find the existence of of a critical mass rn^ ri f)12|) at which the 
persistent current attains a constant value I sa t (USD- This is actually the current produced by 
a massless particle with same charge. The current I sa t appears as a lower cutoff as we further 
decrease the mass. For m < m^ ri the current shows a divergence near the integral values of flux 
(|14p . We also studied the case of a massless fermion. Here also we find that the current has 
the same lower cutoff value (I sa t) and it diverges sharply at the integral flux values. This is a 
possible scenario in graphene ring or on the surface of a topological insulator where the electrons 
behaves as massless Dirac fermions. We then show the nature of total current for a collection 
of rings for different masses. For large mass {figWp) we see the result is similar to that of a 
nonrelativistic particle where total current varies linearly with flux. As we go to the lower mass 
regime the behaviour becomes more nonlinear. The flux-current curve start showing periodic 
plateaus centring at integral flux values {figWfi)- At m = fn^ ri total current becomes a multiple 
step function of flux (fig&jp). The height of the steps are integral multiple of I sa t and each 
individual step is 2I sa t high. This behaviour is solely due to relativistic nature of the particle. 
As we further decrease the mass, the width of the plateau also keeps decreasing (fig\fi$, e). At 
massless limit the total current appears as pulses at integral values of flux (fig^f- One can 
look for such behaviour of total persistent current in a suitable experiment and there occurrence 
would be a direct evidence of the relativistic behaviour of the participating particles. Finally 
we studied the nature of the susceptibility. At large mass the ring behaves as a diamagnet with 
constant negative susceptibility. With decrease of mass it starts showing periodic oscillation 
with negative amplitude and period $o- The ring becomes nonmagnetic (zero susceptibility) 
at critical mass. On further decrease of mass the susceptibility becomes positive. We further 
focus on a particular flux value (<& = 0) and find this variation of susceptibility with mass to 
be continuous. 

In this paper we have discussed everything keeping mass as a parameter. But for experi- 
mental purpose variation of mass is often not a very useful idea. One can alternatively fix the 
mass and vary the radius (R) of the ring. The critical radius (R^) (|15p at which the current 
will become constant (I S at) is found to be half of the Compton wavelength of the particle. But 
in this way we can't go to the massless limit. On the other hand doped graphene can be a good 
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testing ground for such experiment. On a pure graphene sheet electrons behave as massless 
Dirac fermions. In case of a doped graphene there is band gap |31j-|33j. which corresponds a 
massive Dirac fermion. The band-gap can be controlled by the doping concentration. Thus 
instead of mass one can use doping concentration as the driving parameter and study the na- 
ture of the persistent current as well the diamagnet to paramagnet transition for relativistic 
fermions. 
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